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Abstract-Cardiac  reentry,  the  underlying  mechanism 
of  both  tachycardia  and  fibrillation,  is  often  termi¬ 
nated  clinically  with  extracellular  stimuli.  We  devel¬ 
oped  a  one-dimensional  mathematical  model  to  find 
the  probability  that  a  short,  randomly-timed  extracel¬ 
lular  stimulus  would  terminate  reentry  and  to  study 
the  mechanisms  responsible  for  termination.  Our  ring 
model  consisted  of  100  atrial  cells,  with  the  cell  mem¬ 
brane  represented  mathematically  by  a  model  devel¬ 
oped  by  Nygren  et  al.  Stable  reentry  was  estab¬ 
lished,  and  then  a  single  extracellular  stimulus  with 
a  pulsewidth  of  either  2.5  msec  or  5.0  msec  was  ap¬ 
plied  through  a  cathode  positioned  over  cell  25  and 
an  anode  over  cell  50.  Total  delivered  charge  was 
kept  constant.  The  timing  of  the  stimulus  was  varied 
in  1.0  msec  steps  so  as  to  sample  one  complete  rev¬ 
olution  of  the  reentering  wavefront.  The  probability 
that  the  2.5  mec  stimulus  would  terminate  reentry  was 
2.8%,  0.2%  at  the  cathode  and  2.6%  at  the  anode.  The 
probability  of  termination  increased  to  4.2%  with  the 
5.0  msec  stimulus  (1.5%  at  the  cathode  and  2.7%  at 
the  anode).  The  anode  was  more  significant  than  the 
cathode  in  terminating  reentry,  and  the  anodal  mech¬ 
anism  was  less  sensitive  to  the  stimulus  pulsewidth. 
Keywords  -  Electrophysiology,  Reentry,  Electrical 
Stimulation 

I.  Introduction 

Cardiac  reentry  is  believed  to  be  the  primary  mecha¬ 
nism  underlying  the  maintenance  of  tachycardia  and  fib¬ 
rillation  in  the  heart,  both  in  the  ventricle  and  the  atria. 
In  an  attempt  to  terminate  or  modify  reentrant  excitation, 
extracellular  stimuli  are  often  used;  these  interventions 
may  or  may  not  work,  for  reasons  that  are  often  unclear. 
One  difficulty  is  that  a  number  of  seemingly  straightfor¬ 
ward  questions  about  the  interaction  of  the  fields  pro¬ 
duced  by  the  stimuli  with  the  propagating  wavefronts  are 
unanswered  in  the  literature.  For  example,  what  is  the 
probability  that  a  short,  randomly-timed  stimulus  will 
terminate  excitation  in  a  reentrant  loop?  By  what  mech¬ 
anisms  is  termination  possible? 

While  reentry  in  the  whole  heart  is  most  likely  a  three- 
dimensional  phenomenon,  the  simplest  manifestation  is 
continuous  excitation  around  a  thin  ring  of  tissue  (ap- 
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Fig.  1.  Ring  model  consisting  of  100  atrial  cells.  The  diagram  in 
the  center  of  the  ring  shows  the  equivalent  electrical  circuit  for  cell 
0  and  its  two  neighboring  cells  (AC,  atrial  cell;  ES,  extracellular 
space;  Re,  extracellular  resistance;  Rj ,  junctional  resistance;  NM, 
Nygren  membrane).  Extracellular  electrodes  are  positioned  at  cell 
25  and  cell  50  (C,  cathode;  A,  anode). 

proximately  one-dimensional).  Many  investigators  have 
explored  this  type  of  reentry  in  a  wide  range  of  experi¬ 
mental  preparations  [1-5].  These  preparations  can  also 
be  modeled  mathematically,  and,  with  the  aid  of  com¬ 
puters,  the  dynamics  of  a  reentrant  loop  can  be  observed 
in  detail.  Despite  the  advantage  of  providing  for  a  more 
thorough  analysis,  much  of  the  past  theoretical  work  has 
used  simplified  models  of  the  membrane  ionic  currents 
(i.e.  Fitzhugh-Nagumo)  [6-8],  ignored  the  importance 
of  intercellular  coupling  (gap  junctions)  [9,10],  or  relied 
upon  unrealistic  methods  for  stimulation  (i.e.  injecting 
current  intracellularly  as  opposed  to  having  electrodes  po¬ 
sitioned  in  a  restricted  extracellular  space)  [9-13].  While 
much  has  been  learned  about  reentry  and  its  termination 
from  such  models,  it  is  useful  to  compare  those  findings 
to  results  from  more  physiologically  accurate  membrane 
models,  within  a  more  realistic  structure,  and  under  more 
clinically  relevant  stimulation  protocols. 

II.  Methods 

A.  Atrial  Ring  Model 

We  developed  a  one-dimensional  ring  model  consisting 
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Table  I 

Termination  of  reentry  by  an  extracellular  stimulus  for  two  pulsewidths.  The  ranges  of  stimulus  times  that  ter¬ 
minated  reentry  and  the  width  of  the  resulting  termination  windows  are  shown  for  both  cathodal  (rc)  and  anodal 
(7ai,Ta2,Ta3)  stimulation.  Vc  and  Va  are  the  resulting  probabilities  that  a  randomly-timed  stimulus  will  terminate 
reentry  at  the  cathode  and  anode,  respectively.  Vt  =  Vc  +  Va • 


^ stim 

(msec) 

Tc 

(msec) 

Vc 

Tal 

(msec) 

Ta2 

(msec) 

7 ~a3 

(msec) 

Ell 

(msec) 

Va 

Vt 

2.5 

216.8  -  218.4 
(1.6) 

0.2% 

105.2  -  112.6 
(7.4) 

113.0  -  121.4 
(8.4) 

124.0-  126.0 
(2.0) 

17.8 

2.6% 

2.8% 

5.0 

201.1  -  211.4 
(10.3) 

1.5% 

103.8  -  111.7 
(7.9) 

112.1  -  121.2 
(9.1) 

123.3-  125.1 
(1.8) 

18.8 

2.7% 

4.2% 

of  100  atrial  cells  to  study  the  likelihood  that  a  short  ex¬ 
tracellular  stimulus  would  terminate  reentrant  propaga¬ 
tion  (Fig.  1).  The  cell  membrane  was  represented  mathe¬ 
matically  by  a  model  developed  by  Nygren  et  al  [14] .  Each 
cell,  130  pm  long  and  11  /*m  in  diameter,  was  assumed  to 
be  isopotential  [15,16].  The  gap  junctions  connecting  the 
intracellular  space  of  neighboring  cells  were  represented 
by  one  hundred  50  pS  connexons  in  parallel;  this  resulted 
in  a  total  junctional  resistance  Rj  of  200  MCI.  The  total 
extracellular  resistance  Re  between  cells  was  21.1  MCI, 
based  on  an  interstitial  volume  of  0.8  pL  and  an  extracel¬ 
lular  resistivity  of  100  f 1  -  cm.  The  resulting  mathematical 
equations  were  solved  using  an  explicit  method  (Forward 
Euler)  with  a  temporal  discretization  A ts  of  25  psec  ( ts , 
simulation  time). 

B.  Initiating  Reentry 

Reentry  was  initiated  by  stimulating  one  end  (cell  0) 
of  an  100  cell  fiber,  as  the  resulting  wavefront  propagated 
down  the  fiber,  the  two  ends  were  mathematically  con¬ 
nected  to  form  a  ring.  The  wavefront  returned  to  cell 
0,  which  had  recovered  from  its  initial  excitation,  and 
continued  uninterrupted  around  the  ring.  We  allowed 
reentry  to  continue  for  10  seconds  so  that  the  underly¬ 
ing  dynamics  could  reach  a  steady-state.  In  this  time, 
the  wavefront,  propagating  at  1.88  cm / sec,  made  approx¬ 
imately  14.5  revolutions.  The  final  model  state  variables 
(at  ts  =  10.0  sec )  were  recorded  to  use  as  the  initial  con¬ 
ditions  for  subsequent  simulations. 

C.  Stimulation  Protocol 

Simulations,  2.5  msec  in  duration,  were  performed  to 
measure  the  probability  that  a  single  extracellular  stim¬ 
ulus  would  terminate  the  reentrant  wavefront.  The  stim¬ 
ulus  was  applied  through  a  cathode  positioned  over  cell 
25  and  an  anode  positioned  over  cell  50  (Fig.  1).  We 
tested  pulsewidths  Sstim  of  2.5  msec  and  5.0  msec.  The 
stimulus  current  Istim  was  set  to  15400  pA  and  7700  pA 
respectively  to  maintain  a  constant  total  delivered  charge 
Q stim  of  38.5  pC  (Q stim  —  8 stim  *  Rtirn)  •  The  timing  of 
the  stimulus  was  incremented  by  0.1  msec  between  sim¬ 
ulations  until  one  complete  revolution  of  the  reentering 


wavefront  was  sampled.  The  transmembrane  potentials 
from  each  cell  were  recorded  to  determine  whether  the 
stimulus  did  or  did  not  terminate  reentry,  and  at  which 
electrode. 

III.  Results 

Table  I  summarizes  the  results  of  our  simulations.  Both 
cathodal  and  anodal  stimulation  terminated  reentry  for 
each  pulsewidth.  For  a  Sstim  of  2.5  msec,  the  probability 
that  a  randomly-timed  stimulus  would  terminate  reentry 
was  2.8%;  for  a  8stim  of  5.0  msec,  the  probability  in¬ 
creased  to  4.2%.  The  greater  likelihood  that  the  5.0  msec 
stimulus  would  terminate  reentry  was  mainly  due  to  dif¬ 
ferences  at  the  cathode.  However,  for  each  pulsewidth, 
the  probability  of  terminating  reentry  at  the  anode  was 
greater  than  the  probability  of  terminating  reentry  at  the 
cathode. 

A.  Termination  at  the  Cathode 

Fig.  2A  shows  the  timing  of  the  stimulus  relative  to 
the  waveform  at  cell  25  (the  position  of  the  cathode)  that 
resulted  in  termination  of  reentry.  If  the  cathodal  stim¬ 
ulus  was  timed  such  that  the  membrane  in  the  antegrade 
(counterclockwise)  direction  was  refractory  but  the  mem¬ 
brane  in  the  retrograde  direction  was  sufficiently  recov¬ 
ered,  the  result  was  a  single  propagated  secondary  wave- 
front  in  the  retrograde  direction  (unidirectional  block); 
and  when  this  secondary  wavefront  collided  with  and  an¬ 
nihilated  the  primary  reentering  wavefront,  reentry  was 
terminated  and  the  entire  ring  returned  to  a  resting  state. 
Bidirectional  propagation  occurred  if  the  cathodal  stim¬ 
ulus  interacted  with  recovered,  excitable  membrane  in 
both  the  antegrade  and  retrograde  directions;  the  result¬ 
ing  secondary  retrograde  wavefront  annihilated  the  pri¬ 
mary  wavefront  as  above,  but  the  secondary  antegrade 
wavefront  continued  unimpeded  around  the  ring,  essen¬ 
tially  advancing  the  reentry.  If  the  cathodal  stimulus  was 
delivered  over  more  refractory  membrane,  block  occurred 
in  both  directions,  and  the  stimulus  had  no  visible  effect 
on  the  reentering  wavefront. 

The  separate  termination  windows  (rc)  for  each 
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Fig.  2.  (A)  Transmembrane  potential  at  cell  25,  the  location  of  the 
cathode.  The  white  band  extends  to  the  outer  boundaries  of  the 
cathodal  termination  windows  (rc)  for  the  two  stimulus  pulsewidths, 
showing  the  timing  of  the  stimulus  relative  to  the  propagating  wave- 
front  necessary  to  terminate  reentry  at  the  cathode  via  unidirec¬ 
tional  block.  (B)  The  area  in  (A)  outlined  by  the  small  rectangle,  de¬ 
tailing  the  separate  termination  windows  for  the  2.5  msec  stimulus 
(216.8  —  218.4  msec)  and  the  5.0  msec  stimulus  (201.1  —  211.4  msec). 

pulse  width  are  shown  more  clearly  in  Fig.  2B,  provid¬ 
ing  a  graphical  representation  of  the  data  summarized  in 
Table  I.  The  probability  of  termination  at  the  cathode 
with  a  2.5  msec  stimulus  was  significantly  less  than  the 
probability  of  termination  with  a  5.0  msec  stimulus  (0.2% 
versus  1.5%). 

B.  Termination  at  the  Anode 

Fig  3.  shows  the  timing  of  the  stimulus  relative  to  the 
waveform  at  cell  50  (and  the  two  neighboring  cells)  that 
terminated  reentry  at  the  anode,  which  was  positioned 
over  cell  50.  The  highlighted  band  in  Fig.  3  extends  from 
the  beginning  of  the  first  termination  window  (rQi)  for 
the  5.0  msec  stimulus  to  the  end  of  the  last  termination 
window  (r„3)  for  the  2.5  msec  stimulus  (see  Table  I).  If 
the  leading  edge  of  the  wavefront  was  propagating  into 
cell  49,  50,  or  51  at  the  time  the  stimulus  was  delivered, 
the  primary  wavefront  was  blocked  and  reentry  was  ter¬ 
minated;  otherwise,  the  anodal  stimulus  had  little  or  no 
effect  on  reentrant  propagation. 

Unlike  at  the  cathode,  the  probability  of  terminating 
reentry  at  the  anode  was  similar  for  both  pulsewidths 
(2.6%  for  the  2.5  msec  stimulus;  2.7%  for  the  5.0  msec 
stimulus).  The  locations  of  the  three  termination  win¬ 
dows  (ra i,  Ta 2,  Tas)  in  time  for  each  pulsewidth  were  also 
very  similar  (Table  I) .  We  investigated  the  potential  wave¬ 
forms  around  the  anode  at  the  time  of  the  stimulus  in 
more  detail  and  found  that  stimuli  delivered  during  ra i 
prevented  the  excitation  of  cell  49,  stimuli  delivered  dur¬ 
ing  ra 2  prevented  the  excitation  of  cell  50,  and  stimuli 
delivered  during  ra 3  prevented  the  excitation  of  cell  51. 


IV.  Discussion  and  Conclusion 

The  mechanism  by  which  a  properly  timed  depolariz¬ 
ing  (cathodal)  stimulus  can  result  in  the  unidirectional 
block  of  the  elicited  wavefront,  and  the  importance  of 
this  mechanism  in  both  the  initiation  and  termination 
of  reentry,  is  well  represented  in  the  literature  [2,7,8,12, 
13,17-22].  However,  what  is  less  well  understood  is  the 
ability  of  an  anodal  stimulus  to  terminate  reentrant  prop¬ 
agation.  We  found  the  anodal  mechanism  to  be  more 
significant  than  the  cathodal  mechanism,  as  well  as  less 
sensitive  to  the  pulsewidth  of  the  delivered  stimulus. 

These  results  need  to  be  extended  to  include  additional 
stimulus  waveforms.  But  it  is  clear  from  these  limited 
cases  that  developing  protocols  to  more  effectively  ter¬ 
minate  reentrant  propagation  in  the  heart  will  require  a 
more  complete  understanding  of  how  anodal  stimuli  in¬ 
teract  with  the  surrounding  tissue. 
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